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Disconnected diagrams are expected to be sensitive to the inclusion of dynamical fermions. We present a 
feasibility study for the observation of such effects on the nucleonic matrix elements of the axial vector current, 
using SESAM full QCD vacuum configurations with Wilson fermions on 16 3 x 32 lattices, at j3 = 5.6. Starting 
from the standard methods developed by the Kentucky and Tsukuba groups, we investigate the improvement 
from various refinements thereof. 



1. Introduction 

The calculation of hadronic matrix elements 
containing flavor singlet quark bilinear operators 
O = qTq is expected to reveal sea quark effects in 
full QCD simulations. 

The standard method for such calculations uses 
the ratio of the three- and the two-point func- 
tions: 

R(f , W)E n P(n)#(0)) 
H(t) - (N{t)N(0)) 

const + Zo 1 (N\0\N)t, (1) 

where Z$ is the lattice renormalisation factor for 
O. fn the limit of large t R(t) shows a constant 
slope, which determines the matrix element un- 
der investigation. In this work we focus on scalar 
and axial vector insertions (r = 1, T = 7^75) 
needed for the determination of the 7rA^cr term, 
which measures chiral symmetry breaking and for 
the axial coupling constant, g\. In both cases 
connected and disconnected diagrams have to be 
dealt with. For the former suitable methods are 
available. However, the latter present a severe 
bottleneck in form of the computation of A = 
(rM -1 )^^^ which is of complexity ©(volume) [0. 

2. Improving Scalar Insertions 

To tackle the bottleneck several estimator tech- 
niques have been devised by different groups 
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Q. All these methods approximate A by scalar 
products (^rM' 1 ))) with appropriate source 
vectors r\. In the so called wall source method 
a homogeneous 4-volume source is used[|, lead- 
ing to unbiased observables only in the limit of 
an infinite ensemble of gauge configurations. In 
view of the cost of sampling this poses a non triv- 
ial problem for full QCD simulations. 

This problem is alleviated by the stochastic es- 
timator techniques (SET)§,||. 

It has been shown previously |(| that a complex 
Zi noise source is superior to Gaussian noise, in 
the case of the scalar insertion. In this contribu- 
tion we present new methods to improve on the 
signal to noise ratio for scalar and axial vector 
insertions. 

One possible strategy to reduce the noise is to 
restrict the summation over n in eq. g to the 
plateau region of the proton correlator, i.e. to the 
ground state regime of the proton operator. We 
refer to this approach as PSQL (plateau sampling 
of quark loops). 

The efficiency of PSQL (with 1 < n t < t - 1) 
is illustrated in figure [lj where we compare with 
standard SET. Note that the signal to noise ratio 
is improved by a factor 2 — 3. 



3. Axial insertions 

The matrix element of the axial vector current 
between proton states (for zero momentum) can 
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Figure 1. R(t) for the scalar insertion mea- 
sured without PSQL (left) and with PSQL on 200 
SESAM configurations, 100 Z2 noise estimates, 
k = 0.157. 



Figure 2. The standard error cr/y/N on one con- 
figuration for Tr(M ) and Tr ( 7 i 75 M _1 ) urai/i 
Z2 noise source vectors. 



be written as: 



u, d, s. 
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The quantity of interest is the flavor singlet cou- 
pling 
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To determine g\ on the lattice one needs to 
compute Tr(7fc75-M -1 ) which implies off diagonal 
elements of M _1 in spin space [||. Unfortunately 
the SET becomes inefficient for such matrix ele- 
ments, as illustrated in figure ^|. To improve the 
situation we resort to an explicit, component- wise 
treatment in spin space. This amounts to the ap- 
plication of one independent stochastic inversion 
per spin component (spin explicit method, SEM). 
We emphasize that SEM can be combined with 
PSQL. 

The effect is illustrated in figure || and |] for the 
standard error and the mean values on a single 
configuration. We observe earlier asymptotics in 
terms of the number of inversions (the factor 4 
for SEM has been taken into account). 

In order to test the sensitivity of the method 
for a determination of g\ we have analysed 200 
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Figure 3. The standard error a/yN on one 
configuration for Tr^i^M^ 1 ) for Z2 and SEM 
noise source vectors. 



SESAM configurations, at k = 0.157. From the 
ratio 
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Figure 4. Mean values for Tr^i^^M -1 ) on 
one configuration for Z2 and SEM noise source 
vectors. 



4. Conclusion and outlook 

We have presented a variant of the SET which 
is suited to measure the flavor singlet axial vec- 
tor coupling of the nucleon on a limited sample of 
QCD configurations. We can cope with the small- 
ness of the axial vector signal by spin explicit in- 
versions. Although this amounts to an increase 
in computational cost by a factor 4 for each esti- 
mate, the overhead turns out to be largely com- 
pensated by earlier asymptotics. The second in- 
gredient of our proposal is the use of plateau sam- 
pling, which appears to reduce the fluctuations by 
a factor 2 — 3, both for the scalar and for the axial 
insertions. 

This encourages us to perform a full physics 
analysis to the complete set of SESAM configu- 
rations. 



we obtain signals as given in figure |. Using 
both, SEM and PSQL, one finds the signal to 
stick clearly out of the noise! 
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Figure 5. R(t)^ c measured without PSQL (left) 
and with PSQL on 200 SESAM configurations, 
100 SEM noise source vectors, n — 0.157. 



